Symmetry and disorder of the vitreous vortex lattice in an overdoped BaFe 2 _ x Co x As 2 
superconductor: Indication for strong single-vortex pinning 
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The disordered flux line lattice in single crystals of the slightly overdoped BaFe 2 _ x Co ;c As 2 (x = 0.19, 
T c = 23 K) superconductor is studied by magnetization measurements, small-angle neutron scattering 
(SANS), and magnetic force microscopy (MFM). In the whole range of magnetic fields up to 9T, vortex 
pinning precludes the formation of an ordered Abrikosov lattice. Instead, a vitreous vortex phase (vortex 
glass) with a short-range hexagonal order is observed. Statistical processing of MFM datasets lets us directly 
measure its radial and angular distribution functions and extract the radial correlation length £. In contrast 
to predictions of the collective pinning model, no increase in the correlated volume with the applied field 
is observed. Instead, we find that £ decreases as (1.3 ± 0.1)^ oc H~ 1/2 over four decades of the applied 
magnetic field, where Rj is the radius of the first coordination shell of the vortex lattice. Such universal 
scaling of £ implies that the vortex pinning in iron arsenides remains strong even in the absence of static 
magnetism. This result is consistent with all the real- and reciprocal-space vortex-lattice measurements in 
overdoped as-grown BaFe 2 _ x Co, t As 2 published to date and is thus sample-independent. The failure of the 
collective pinning model suggests that the vortices remain in the single-vortex pinning limit even in high 
magnetic fields up to 9 T. 

PACS numbers: 74.70.Xa 74.25.Wx 68.37.Rt 61.05.fg 



I. INTRODUCTION 

The behavior of supercurrent vortices associated with 
magnetic flux quanta penetrating a type-II superconductor, 
known as Abrikosov vortices, 1 contains valuable information 
about the microscopic material properties, such as the sym- 
metry of the superconducting order parameter or the vortex 
pinning force. 2 ' 3 These, in their turn, determine essential 
macroscopic physical characteristics of the superconductor, 
such as the critical current density or critical fields. In clean, 
highly homogeneous superconductors the vortices form an 
ordered Abrikosov lattice due to the repulsive intervortex 
interactions, as was first observed experimentally by neu- 
tron diffraction in 1964 4 and later imaged directly by Bitter 
decoration. 5 ' 6 At low fields, the close-packed triangular flux- 
line lattice is commonly realized. As the field is increased, 
however, the strong anisotropy of the Fermi surface 7 and of 
the superconducting gap 8 in some unconventional supercon- 
ductors may lead to pronounced phase transitions to less 
trivial lattice symmetries as soon as the anisotropic vortex 
cores are forced into closer proximity. 9,10 ' 11 ' 12 ' 13 

In many actual materials, nevertheless, pinning of the 
magnetic flux lines by point-like defects, disorder, or 
twin domain boundaries prevents the formation of a well- 
ordered vortex lattice. Instead, a vortex-glass phase char- 
acterized only by a short-range positional order is often 
formed. 2 ' 14 ' 15 ' 16 ' 17 The presence of such disorder is not nec- 
essarily an issue of sample quality. High-T c superconductors, 
for example, are obtained by chemical doping of stoichio- 



metric parent compounds with impurity atoms or vacancies, 
which by themselves may serve as natural pinning centers 18 
and enhance pinning even in pristine single crystals. 

There are two complementary classes of experimental 
techniques for visualization of magnetic-flux-line lattices. 
The first one includes diffraction methods like small-angle 
neutron scattering (SANS), capable of observing the vor- 
tex lattice in the reciprocal space by measuring neutron 
diffraction patterns that originate from the inhomogeneous 
field distribution in the whole volume of the sample. The 
second class includes real-space imaging methods, 19 such as 
Bitter decoration, scanning tunneling microscopy (STM), or 
magnetic force microscopy (MFM), that can directly probe 
the local distribution of vortices within a small area on the 
sample surface. 

In the present paper, we apply both SANS and MFM as 
complementary techniques to study the magnetic-flux-line 
lattices in slightly overdoped BaFe 2 _ x Co x As 2 (BFCA) with 
the doping level x = 0.19 and the onset superconducting 
transition temperature T c = 23 K. This material belongs to 
the so-called 122-family of the novel iron arsenide high- 
temperature superconductors. 20 ' 21 ' 22 Observations of vortex 
lattices in BFCA are still limited to a recent STM study on an 
overdoped (x = 0.2) sample, 23 a combination of SANS and 
Bitter decoration methods applied to a slightly underdoped 
(x = 0.14), 24 ' 25 and an MFM investigation of an underdoped 
(x = 0.10) 26 single crystal. All these studies revealed a 
highly disordered flux-line lattice that is suggestive of strong 
bulk pinning, in agreement with our results presented below. 



In the following, we will quantify the degree of disorder 
in such a lattice using statistical analysis of multiple MFM 
images. We will show that the vortex lattice is characterized 
by predominantly hexagonal coordination, as expected for 
a disordered triangular lattice, and will estimate the radial 
correlation length of such order from our MFM and SANS 
data in a wide range of the applied magnetic fields. 



n. MODEL OF THE LATTICE DISORDER 

Comparing the correlation length of the vortex lattice 
extracted from SANS data, measured in the reciprocal space, 
to the direct-space MFM measurements is rather nontrivial, 
as different models of disorder may lead to ambiguous defi- 
nitions of the correlation length. Strictly speaking, no single 
well-defined parameter can fully characterize the "degree of 
disorder" in vitreous or amorphous lattices with short-range 
correlations. Therefore, to gain some consistent quantitative 
understanding from both direct- and reciprocal-space mea- 
surements, we start by introducing a continuous Gaussian 
model of positional disorder applied to flux line displace- 
ments perpendicular to the field direction. This will later 
prove to be consistent with our experimental results. 

Let r ;j be a set of vectors representing a regular periodic 
two-dimensional (2D) lattice in direct space with a lattice 
parameter a. The disorder is introduced by adding a random 
vector w(r i; ) = {w x (r ;j ), w y (r ; ,)} to every lattice site: 



Ry = ry+w(r y ). 



(1) 



The vector components w x (r) and w y (r) are two inde- 
pendent centered Gaussian stochastic processes of a two- 
dimensional argument with dispersion given by 

Vr X) r 2 eR 2 ([w^rj - w x (r 2 )] 2 } 

= ( KCrJ - w y (r 2 )] 2 ) = a 2 \r t - r 2 |/a. (2) 

Here o <C a is the standard deviation of the distance be- 
tween nearest-neighbor lattice sites from its mean value a. 
Angle brackets denote mathematical expectation. 

The radial distribution function (RDF) 27 ' 28 of such a lattice, 
g(r), can be well approximated by a sum of Gaussian peaks, 
such that their widths increase proportionally to Jr, as 
follows from the property (2) : 
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Here R n is the radius of the n th coordination shell, and N n is 
the number of sites in this shell. 

Under these assumptions, on the scale of L = a 3 /(2a) 2 
the width of the peaks, 2a^/ r/a, will become comparable 
with the lattice spacing, and individual peaks in g(r) will 
merge into a linearly increasing continuum background. 
The knowledge of L therefore provides an upper boundary 
for the range of order — a distance beyond which no mea- 
surable oscillations in g(r) can be detected. 27 For r < L, 
the RDF will exhibit oscillatory behavior, with the roughly 



exponentially decaying amplitude ~ exp(— r/Q, as we will 
demonstrate below. Such exponential decay suggests a nat- 
ural definition for the radial correlation length (. We point 
out that £ is several times smaller than the maximal length 
scale at which detectable correlations exist. More precisely, 
the following relationship between £ and the parameters of 
our model is derived in the Appendix: 
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The reduced radial distribution function g(r) can be con- 
structed as 29 



gM = g(r) - 2urg , 



(5) 



where g is the average density of sites in the 2D lattice. For 
a triangular lattice, in particular, g = 2/\/3a 2 . The new 
function g(r) has a convenient property 72 : 



lim g(r) = 0, 

r— »oo 

so that its Fourier transform is well defined: 



S(fc): 



g(r)sin(fcr)dr. 



(6) 



(7) 



It can be shown that S(fc) is nothing else but the reduced 
structure factor, which can be obtained from the experimen- 
tally measured elastic scattering cross-section after form- 
factor and resolution corrections. 29 - 30 - 31 The exponential 
decay of oscillations in g(r) ensures a Lorentzian shape 
of the first diffraction peak in S(fc) with the half- width at 
half-maximum (HWHM) of £ _1 (for example, see Ref. 32). 

IE. SAMPLE PREPARATION AND CHARACTERIZATION 

The single crystal of BFCA used for the present study has 
been grown by the self-flux method in a zirconia crucible, us- 
ing an A1 2 3 nucleation center. 33 SANS measurements were 
done on a ~ 0.4 g as-grown single crystal, whereas a smaller 
sample from the same growth batch was used for MFM 
measurements. The doping level of x = 0.19 ± 0.01 was de- 
termined by energy-dispersive X-ray analysis (EDX) after av- 
eraging over the sample surface. The onset superconducting 
transition temperature T c = 23 K was measured by SQUID 
magnetometry [Fig. 1 (a)]. These values let us pinpoint our 
sample on the overdoped side of the phase diagram. 34,35 ' 36 
The choice of such doping level is motivated by the non- 
magnetic ground state of this system, which excludes the 
possibility of any phase separation 37 ' 38 ' 39,40 ' 41 or microscopic 
coexistence 42,43,44 of the antiferromagnetic and supercon- 
ducting phases that is observed in underdoped samples and 
can potentially influence the flux penetration pattern. Since 
it was proposed that in the underdoped single crystals vortex 
pinning occurs mainly on the boundaries of intertwinned or- 
thorhombic/ antiferromagnetic domains, 45,46 ' 47 by choosing 
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Fig. 1 (color online). Magnetization measurements, (a) dc-susceptibility measured after zero-field cooling in various applied fields 
up to 7 T. (b) dc susceptibility measured during field-cooling and -warming. In both panels, note the different scale for the 1 mT and 
lOmT data (right axis) and the IT, 3T, and 7T data (left axis), (c) Isothermal magnetization loops measured at various temperatures (as 
indicated beside each curve) in a wide field range up to 14 T exhibiting the "fishtail" effect. Note that the position of the second peak H p 
is somewhat higher in measurements with increasing field than in the reverse direction, (d) Magnetic phase diagram summarizing the 
temperature dependence of the second magnetization peak H p , irreversibility field H irr , and second critical field H c2 extracted from the 
magnetization measurements. The values of H p measured in increasing and decreasing fields are marked with up- and down-pointing 
triangles, respectively. Triangles and squares denote values of H in and H c2 measured on two different sample pieces from the same growth 
batch. 



an overdoped sample we can exclude this pinning channel 
from consideration. The absence of static magnetism in our 
sample has been confirmed by a zero-field low-temperature 
muon spin relaxation (/xSR) experiment. 

Temperature-dependent magnetization measured by 
SQUID magnetometry after zero-field cooling (ZFC) and 
during in-field cooling (FC) with subsequent warming are 
shown in Fig. 1 (a) and (b), respectively. The temperatures 
at which the FC curves start to deviate from linear behav- 
ior correspond to the superconducting transition. The dip 
at which the hysteresis develops represents the irreversibil- 
ity temperature T irr , below which the flux line lattice gets 
pinned, whereas at temperatures between T irr and T c de- 
pinning of the vortices by thermal fluctuations leads to the 
melting of the flux line lattice. Combining the ZFC/FC data 



for various applied fields, the temperature dependencies of 
the second critical field H c2 and of the irreversibility field 
H irr , were extracted. 

Isothermal magnetic hysteresis loops were measured at 
different temperatures in a wide range of magnetic fields 
up to 14 T using a Quantum Design PPMS extraction sample 
magnetometer, as shown in panel (c) of the same figure. 
The non-monotonic behavior with a well-pronounced sec- 
ond peak in magnetization represents the so-called "fish- 
tail" effect. It is well known from both conventional and 
high-temperature superconductors, but its origins are dis- 
cussed controversially. Most theoretical approaches agree 
that the temperature-dependent field H p , at which the sec- 
ond magnetization peak has its maximum, is related to 
vortex pinning and corresponds to a crossover between 
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two different regimes of the vortex lattice. Various scenar- 
ios have been proposed responsible for the fishtail effect, 
among them a change in dynamics of the vortex lattice, 48 
a change in flux creep behavior, 49 or a change in its elas- 
tic properties. 50 A fishtail effect of comparable magnitude 
has already been observed recently in several iron arsenide 
superconductors. 46,47 ' 51 ' 52 - 53 

The values of all three characteristic fields (H c2 , H irr , and 
Hp) define the magnetic phase diagram of the material that 
is summarized in Fig. 1 (d). The superconducting phase 
consists of three regions: (i) vortex liquid phase above H irr , 
where magnetic flux lines are not pinned due to strong ther- 
mal fluctuations; (ii) region between H p and H irr , where 
intervortex repulsion exceeds the typical pinning forces, pre- 
sumably leading to an elastically interacting vortex lattice; 
(iii) low-field region, the nature of which may depend on the 
actual pinning strength. 50,54 ' 55 The measurements presented 
in the following were performed at low temperatures in a 
broad range of magnetic fields spanning the region below 
Hp to elucidate the actual structure of the vortex lattice 
within this phase. 

IV SMALL-ANGLE NEUTRON SCATTERING 

A. Experimental conditions 

The SANS experiment was carried out using the NG3- 
SANS diffractometer at NIST. For all measurements pre- 
sented here, cold neutrons were collimated over a distance 
of 8.0 m before the sample. The diffracted neutrons were col- 
lected using a two-dimensional multidetector of 128 x 128 
pixels, each 5.08 mm x 5.08 mm, positioned 11.5 m behind 
the sample. 

The momentum resolution of the instrument Aq(q), ex- 
pressed in terms of the Gaussian standard deviation, was 
estimated in the approximation of strong longitudinal disor- 
der, i.e. large rocking curve width. In this limit, the instru- 
mental contribution to the peak width in the radial direction 
is given by 56 

Aq(q)= >/o- a 2 fc 2 + q 2 (Afc/fc) 2 (8) 

For our SANS results herein, the standard deviation angu- 
lar spread of the incoming beam cr a was measured to be 
0.085°, the mean neutron wave number fc was 1.26A -1 
with standard deviation spread Afc/fc = 6.2 %. 

The single-crystalline sample was loaded into the cold 
variable-temperature insert of a cryomagnet that provided 
applied fields of up to 9 T and a base temperature of 3.5 K. 
The sample was oriented with the a-axis vertical, i.e. per- 
pendicular to the applied field, and the c-axis parallel to 
the field. Our measurements were carried out in the usual 
SANS experimental geometry, with applied field approxi- 
mately parallel to the neutron beam. The vortex lattice 
was prepared in the sample by applying the desired field 
above T c , and subsequent field-cooling to 3.5 K. In all cases, 
background measurements were carried out at 26 K (above 
T c ) and subtracted from the field-cooled foreground mea- 
surements. An attempt to oscillate the field during cooling 



to equilibrate the vortex lattice did not result in any notable 
reduction of lattice disorder. Since such wiggling helps only 
against a weak pinning background, this fact is consistent 
with strong vortex pinning in our sample. 

B. SANS results 

Usually in SANS experiments, the sample and cryomagnet 
are rotated together in order to allow a reciprocal lattice 
vector of the vortex lattice to satisfy the Bragg condition 
at the detector. Generally speaking, the Bragg spots in the 
reciprocal space have a finite size determined both by the 
instrument resolution and by the structure of the vortex 
lattice. Therefore, in order to capture all of the diffracted 
intensity associated with a Bragg spot, the sample is usually 
rotated through a range of angles about the vertical axis 
(the rocking scan) . However, instead of observing distinct 
Bragg spots, our preliminary measurements on the BFCA 
sample identified the vortex lattice as strongly disordered; 
the diffraction signal was weak and took the form of a ring 
centered on the origin of reciprocal space. This observation 
is qualitatively identical to the previous studies reported by 
Eskildsen et. al. on a slightly underdoped sample. 24 ' 25 On 
rotating the sample by ±6° about the vertical axis, within 
the statistical error there was no variation of the diffracted 
intensity at any point around the ring. This indicates that 
the rocking curve width is extremely broad, and hence all 
subsequent measurements were carried out at a fixed rota- 
tion angle. The observation of a ring of scattering indicates 
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Fig. 2 (color online) . Angle-averaged SANS intensity as a function 
of momentum transfer |q|, measured at different magnetic fields 
between 0.9 T and 2.0 T. For clarity, each curve is offset by an 
increment of 0.3 units from the one below it. The vertical scale 
represents the mean diffracted intensity per pixel, normalized to 
the standard monitor. The solid lines are Gaussian fits to Eq. (9) . 
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Table I. Best-fit parameters of the SANS data from Fig. 2. Here a A = y 2(p /H</3 and a n = sj <p /H stand for the parameters of perfect 
triangular and square vortex lattices, respectively, calculated for every magnetic field, <j> = h/2e being the magnetic flux quantum; 
4Ti/a A ^/3 and 2n/a n are the expected radii of the corresponding diffraction rings that are to be compared with the fitted peak position 
q ; HWHM of the diffraction peaks are given prior to resolution correction; ( is the resolution-corrected radial correlation length. 



the absence of orientational order of the vortex lattice about 
the field axis. Both of these observations are consistent 
with strong bulk pinning of the vortices. However, the finite 
radial width of the diffraction ring indicates a finite planar, 
or d -spacing, order which we can investigate in a more 
quantitative manner. 

In Fig. 2, for a selection of fields we show the radial 
dependence of the diffracted intensity. At each field, the 
data were obtained by averaging over 2tc annuli that were 
centered at the origin of reciprocal space, as defined from 
the position of the unscattered beam on the detector. For 
each |q|-bin of a curve, the intensity is obtained from the 
mean counts per standard monitor per pixel averaged over 
all the pixels whose location in q lies within the width of the 
|q|-bin from the origin. The resulting curves show weak, but 
clear, peaks in the radial intensity, where the location of the 
peak in |q| corresponds to the radius of the diffraction ring 
in reciprocal space. With increasing field, we clearly see the 
position of the peak in |q| increasing, thus confirming that 
the origin of our diffraction signal is due to a vortex lattice. 
The solid lines shown in Fig. 2 are fits to the following 
Gaussian function: 

/(q)=Aexp[-(q-q ) 2 /2(r 2 + Aq(q ) 2 )] . (9) 

The fitting results are summarized in Table I. 

Despite our relaxed resolution setup, the observed diffrac- 
tion peaks are not resolution-limited. Because of these broad 
widths, it is difficult to draw a confident conclusion about 
the local symmetry of the vortex lattice solely from the SANS 
data. Within the uncertainty of our results, the radius of the 
diffraction ring could be reconciled both with the hexagonal 
and square coordination (see Table I). 

In principle, the broad peaks could indicate that the vor- 
tex lattice is composed of very small and randomly ori- 
ented anisotropic domains. Anisotropic vortex structures, 
for example in the form of distorted hexagonal vortex lat- 
tices, are commonly observed in SANS studies of type-II 
superconductors. 9,10,11,12 ' 13 They are particularly prominent 
at higher fields when the vortices start to overlap. In these 
situations, anisotropies of the superconducting gap, 8 or of 
the Fermi surface in combination with non-local effects, 7 
can directly cause the deviation from the ideal isotropic 
and hexagonal vortex lattice. However, any influence of a 
superconducting gap anisotropy is not expected to be impor- 
tant over the field range explored here. 8 Furthermore, the 



anisotropy caused by non-local effects is expected to be sup- 
pressed in the presence of disorder, due to a corresponding 
reduction in the mean free path (reduction in non-locality 
range). As all the available direct-space measurements, 
both in previous works 23 ' 24,25 and this study, point towards 
isotropic and hexagonal vortex lattice coordinations, we 
assume the lattice disorder to be the dominant source of 
diffraction-line broadening in our SANS data. Therefore, 
the HWHM of the Gaussian lineshapes shown in Fig. 2 can 
be used to measure the radial correlation length £ according 
to Eq. (9), ignoring possible structural contributions. Sub- 
sequent comparison of the extracted values of £ with those 
extracted from direct-space STS measurements performed 
in somewhat higher fields 23 will show a reasonable agree- 
ment between the two probes, indicating the validity of this 
approximation. 

V MAGNETIC FORCE MICROSCOPY 
A. Experimental conditions 

Low-temperature MFM measurements were conducted us- 
ing an Omicron Cryogenic SFM scanning force microscope 57 
supplied with a commercial SSS-MFMR magnetic tip from 
NanoAndMore GmbH with a force constant of ~ 2.8 N/m and 
a resonance frequency of 51 kHz. The measured MFM signal 
represents a shift of the cantilever's resonance frequency, 
which is proportional to the normal component of the force 
gradient acting between the tip and the sample at a given 
scanning distance above the surface. 58 ' 59 

A piece of the same single-crystalline sample as the one 
used for SANS measurements was cleaved in situ in high 
vacuum to get a clean surface. Magnetic fields H = 3 and 
6 mT were applied along the crystallographic c-axis of the 
sample in the direction antiparallel to the magnetization 
of the MFM tip. In such geometry, the interaction between 
the MFM tip and the superconducting vortices is repulsive 
and shows up as a positive frequency shift. The sample was 
then field-cooled down to the base temperature of 8 K. All 
measurements were done at a distance of 80 nm between 
the tip and the sample. Due to a technical limitation of the 
microscope, the size of a single MFM image could not be 
made larger than 4/xm x 4/xm at the base measurement 
temperature. Therefore, to gather enough statistics, 8 and 4 
images were measured at 3 and 6 mT, respectively, with full 
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Fig. 3 (color online) . (a) One out of eight MFM images measured at 3 mT with vortex positions (red dots) determined by 2D Gaussian 
fitting. The lighter colors at the vortex positions correspond to a positive change of the cantilever's resonance frequency, (b) Delaunay 
triangulation of the vortex lattice (black lines) and the corresponding Voronoi tessellation (gray lines) . The vortices are marked according 
to the number of nearest neighbors, (c) Coordination number distribution for all 3 mT images, showing the predominance of a hexagonal 
coordination. The thin line is a Gaussian fit of the histogram. 



temperature cycling before each measurement. 

B. Statistical processing of MFM datasets 

For each vortex, the 2D vector of its coordinates v, was 
determined by 2D Gaussian fitting, the width being kept 
equal for all the peaks. This method allows for sub-pixel 
precision and can successfully locate even those vortices 
that are partially cut by the image boundary, as shown in 
Fig. 3 (a) . Voronoi cells have been calculated for all vortices, 
i.e. the loci of points that are closer to v, than to any other 
vortex in the set. 60 ' 51 The number of sides of each Voronoi 
cell therefore represents the number of nearest neighbors 
for the corresponding vortex. Straight lines connecting all 
nearest-neighbor sites form a Delaunay triangulation, 62 as 
illustrated in Fig. 3 (b). Panel (c) of the same figure shows 
the distribution of vortex coordination numbers for all im- 
ages measured at 3 mT. Note that only vortices surrounded 
by closed Voronoi polygons possess the full set of neighbors 
and can therefore be counted. The maximum of the distri- 
bution is reached for the coordination number of 6, which 
gives us the first indication of the local lattice symmetry. 
In the following, we will therefore assume that the vortex 
phase can be treated as a disordered triangular lattice. More 
solid evidence supporting this assumption will come from 
the consideration of the angular distribution function (ADF) 
in what follows. 

In the next step, the relative distances dy = |v f — V; | were 
determined for all pairs of vortices v, and v,-. The histograms 
of such distances, which approximate the RDF, are plotted 
in Fig. 4 (a) for both magnetic field values. Solid lines show 
fits to 

ijVg(r)5r 
l + exp[(r-R max )/W max ]' 

where N is the total number of vortices in the dataset, g(r) 
is the RDF of a triangular lattice given by Eq. (3), 5r is 
the histogram's bin size, and the factor of 1/2 is added 



to avoid double counting of equivalent vortex pairs (dy = 
dji). The denominator is an empirical factor introduced to 
compensate for the distribution function cut-off at large r 
due to a finite image size. The best-fit parameters of the 
model are given in Table II. The pristine RDFs deprived 
of the experimental correction factors are also plotted in 
Fig. 5 (a). 

The most pronounced peak appearing at r = R x corre- 
sponds to the first coordination shell of the vortex lattice and 
can be clearly seen in both histograms. A good agreement is 
found between the radius of the first coordinatio n shell and 
the expected lattice constant a = \J 2<p /HV3 calculated 
for a perfect triangular vortex arrangement (see Table II). 
Here <£ = h/2e = 2.07 x 10" 15 Tm 2 is the magnetic flux 
quantum. The ratio L = Rj/(2cr) 2 that provides an upper 
boundary for the range of order in the lattice is the same 
for both fields within the statistical error. The direct esti- 
mation of the correlation length from the RDF, illustrated 
in Fig. 5 (b), also results in the value of £ = 1.1 ± 0.2 /xm 
for both fields. This value was obtained by extracting the 
reduced RDF from the experimental fits, plotting its absolute 
value in logarithmic coordinates and fitting the positions of 
its numerous maxima with an exponential-decay function 
expfA — r/Q, as shown by the dashed lines in the figure. 

Note that the resulting value of £ is nearly a factor of 5 
smaller than the range of order L, in agreement with Eq. (4). 
Therefore, distinct correlation peaks can be still seen in the 
RDF histograms [Fig. 4 (a)] even at distances of about 2 to 
3 correlation lengths, depending on statistics. For example, 
the second peak appearing around 1.6 /im in the 3 mT data 
and around 1.2 /jm in the 6mT data corresponds to the un- 
resolved second and third coordination shells of a triangular 
lattice, as indicated by the corresponding coordination shell 
radii shown by vertical lines. The height of these lines is 
proportional to the number of sites in each shell. 

Next, we study the angular distribution of the vortex 
positions. For this, all possible unordered vortex triples 
{v,, v,-, VjJ are considered in every MFM image, such that v.- 
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H = 3mT H = 6 mT 




Fig. 4 (color online) . Statistical analysis of multiple 4 jim x 4 jum MFM images measured in magnetic fields of 3 mT (left) and 6 mT 
(right). For the two fields, data from 8 and 4 images, respectively, were used, (a) Radial distribution function. Vertical lines mark the 
positions and relative intensities of individual coordination shell peaks, (b) Triplet distribution function, (c) Angular distribution function. 
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H 


N 


a A = ^24> /HV3 


■^1 ® ^max ^rnax 


C L=R\/{2a? 


(deg) 


3.0 mT 
6.0 mT 


166 
163 


0.893 jum 
0.631 [im 


0.90 0.19 1.72 0.81 
0.65 0.11 1.82 0.73 


1.1 ±0.2 5.1 ±0.4 
1.1 ±0.2 5.7±0.4 


15.1 ± 1.0 

14.2 ± 1.0 



Table II. Best-fit parameters of the RDF and ADF extracted from MFM data and the derived values of the radial correlation length £ and 
the range of order L. All values are given in /zm, unless otherwise specified. The lattice constant a A , shown in the second column, is the 
expected value calculated for a perfect triangular lattice. It is to be compared with the measured radius of the first coordination shell R 1 . 
The statistical errors of £, L, and cr^ were estimated as the standard deviation of values obtained for different binning of the histograms. 



lies within the first coordination shell from V; . The threshold 
for the distance between vortices to achieve this limitation 
was chosen to include most of the first-shell peak in RDF 
(i.e. somewhat above the maximum of the peak). The 
optimal threshold value for our case was found by trial 
and error to be ~ 1.0 /im at 3 mT and ~ 0.75 /im at 6 mT. 
For all such triples, we plot a point in polar coordinates 
[Fig. 4(b)] at the distance |v fc | from the origin and at an 
angle (p formed by the vectors v.- — v ; and — v ; . Each point 
is plotted twice: for 4> and —cp. The dashed circle shows 



(b) 




H = 3mT,£ = 1.lMm 

H = 6mT, ; = 1.1pm 




Fig. 5 (color online). Determination of the radial correlation 
length £ from MFM data, (a) Radial distribution functions g(r) 
extracted from the fits of the histograms shown in Fig. 3 (a) . (b) Ab- 
solute value of the reduced RDF estimated as 

g(r) = g(r) - 4ur/v / 3Ri - lim [g(r) - 4nr/v / 3Ri] , 

plotted in logarithmic scale to emphasize the exponential decay of 
the oscillations. The slope of the exponential envelopes (dashed 
lines) is a measure of the radial correlation length £, which in our 
case is the same for both fields within the statistical error. 



the threshold of the first coordination shell. In such angle- 
resolved representation, the density of points represents the 
triple-distribution function (TDF), so that diffuse clusters 
of points forming around regular lattice sites can be seen. 
The above-discussed RDF g(r) can essentially be obtained 
by angle-integration of the TDF. 

Now it is straightforward to obtain the ADF of the first 
coordination shell by plotting a histogram of those angles 
4> that correspond to points lying within the threshold ra- 
dius from the origin (i.e. within the dashed circle). Such 
histograms are shown in Fig. 4 (c). They immediately reveal 
a hexagonal coordination of the lattice, as indicated by the 
pronounced peak located at ±60° followed by a dip at ±90°. 
This lets us conclude that the vortex glass phase can be con- 
sidered as a highly disordered triangular lattice. By analogy 
with Eq. (3), the function used to fit the ADF histogram is a 
sum of Gaussian peaks of the form 

q W = ^V - 1 (11) 
^V2n ^ y/2sm(nn/6) 

In this model, cr^ y/2a/a is the standard deviation of 
the first peak located at ±60°. The dispersion of the two 
other peaks at ±120° and 180° is assumed to increase as 
2(7? sin(<£ /2), i.e. proportionally to the distance between 
endpoints of two unit vectors separated by the angle 4>, in 
accordance with Eq. (2) . The correction coefficient A < 1 is 
introduced to compensate for the finite image size effects, 
and 50 is the histogram's bin size. 

The only physically meaningful fitting parameter in this 
model is a^, being just another measure for the degree of 
randomness in the lattice. The convenience of this parame- 
ter is that it is dimensionless, unlike a, £, or L, so it does 
not require comparison to other characteristic length scales 
in the system. Again, the values of cr^ extracted from the 
ADF histograms (see Table II) turn out to be nearly the same 
for 3 mT and 6 mT within the statistical error. 

Here we note in passing that the presented method of 
ADF analysis is very efficient for discerning the local symme- 
try of a short-range order from direct-space measurements. 
Because the position of the first peak in such a plot differs by 
as much as 30° for the triangular and square lattices, which 
in our case corresponds to ~ 2 o^, it provides the optimal 
way to distinguish between the two symmetries. For compar- 
ison, the ratio of the lattice constants for the triangular and 

square lattice symmetries is just a A /a n = y 2 » 1.075. 
Detecting the difference of 7.5% in the radius of the first 



8 



E 

3. 



0.1 



0.01 



BD a » 

r=6K 

MFM 
T=8K 



» BaFe 181 Co 019 As 2 [this work] 
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Fig. 6 (color online) . Low-temperature radial correlation length 
£ of the vortex glass phase summarized as a function of the mag- 
netic field H. Values extracted from the high-field STS data of 
Y. Yin et al. 23 and from the SANS and Bitter decoration (BD) data 
of M. R. Eskildsen et al. 24,25 are added for comparison. Within the 
shaded region the vortex-lattice constant a becomes comparable 
with the penetration depth 2X. The width of this region reflects an 
uncertainty of the reported penetration depth values. 63,64,65 A fit 
to the power law £ oc H~ 1/2 is shown by the dashed line. 



coordination shell R x would be impossible under the same 
conditions. 

VI. SUMMARY AND DISCUSSION 

The radial correlation length £ of the vortex glass phase 
is plotted in Fig. 6 as a function of the applied magnetic 
field ranging from 1 mT to 9 T. In addition to the values 
extracted from our MFM and SANS data, the figure also 
includes correlation lengths derived from the high-field STS 
data of Y. Yin et al. 23 as well as the SANS and low-field Bitter 
decoration data of M. R. Eskildsen et al. 24 - 25 STS and Bitter 
decoration data were analyzed by the same above-described 
algorithm that we used for statistical MFM image processing. 
All presented data were measured on samples with doping 
levels similar to ours and reveal a general power-law trend 
in the field-dependent correlation length that persists over 
nearly four decades of the applied magnetic field: 



£ w (1.3 ± 0.1)7?! = (1.3 ± 0.1)a a H 



-1/2 



(12) 



Such universality implies that the disorder of the flux line 
lattice is sample-independent and can possibly be considered 
an intrinsic property of the as-grown BFCA superconductors 
close to the optimal doping. The good agreement between 
correlation length values measured by direct- and reciprocal- 
space probes confirms the consistency of the two methods. 

Our observations are difficult to reconcile with the col- 
lective pinning model, 66 according to which an increase of 
the correlated volume with the applied magnetic field is 
expected 48,67 that would lead to deviations from the simple 



power law scaling of £ at high fields. Since the correlation 
radius R c (transversal size of the vortex bundle) predicted in 
the collective pinning model is inversely proportional to the 
mean square of the pinning force, 67 in the limit of strong 
pinning it should become smaller than the vortex-lattice 
constant, leading to the breakdown of the collective model. 
In this limiting case of single-vortex pinning, vortices no 
longer form bundles, but behave individually. The ratio £/a 
then remains at a saturated value of the order of unity, in 
agreement Fig. 6. Thus, the observed field dependence of 
the radial correlation length can be understood in the limit 
of strong single-vortex pinning, 50 ' 54 meaning that the typical 
pinning forces dominate over the intervortex repulsion 50 in 
the studied field range. Persistence of such behavior up to 
the fields comparable with H p at low temperatures suggests 
that the second maximum in magnetization possibly repre- 
sents a crossover from the single-vortex pinning regime to 
an elastically interacting vortex lattice. 55,68 

For comparison, in the low-T c superconductor SnMo 6 S 8 , 
where the vortex glass phase has been imaged recently in di- 
rect space by STM, 69 the temperature-dependent crossover 
from the single to collective pinning regimes was observed 
at much lower magnetic fields below H X (T —> 0) = 2.5 mT, 
well separated from H p . In addition, the positional correla- 
tion functions in this material exhibited a slower power-law 
decay, expected for a weakly disordered Bragg glass, 19,70 in 
contrast to the exponential decay that we observe in BFCA. 

The origin of such anomalously strong pinning in doped 
iron arsenide superconductors and the ways to influence 
it artificially remain to be investigated. In a recent work, 
Prozorov et aZ. 46 argue that in an underdoped BFCA the 
pinning occurs mainly on the boundaries of intertwinned or- 
thorhombic/antiferromagnetic domains. 45 Our work, how- 
ever, indicates that pinning still remains anomalously strong 
even in the overdoped single crystals, where static mag- 
netism is fully suppressed. On the one hand, this could 
be an evidence in favor of a purely chemical origin of the 
pinning landscape. On the other hand, evidence for some 
kind of magnetic order induced by weak magnetic fields 
in overdoped iron arsenides has been reported from ,uSR 
experiments. 71 Therefore the influence of magnetism on the 
pinning potential cannot be fully excluded. 
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Appendix 

To derive the relationship (4) between the parameters of 
our disorder model and the radial correlation length £, let 
us consider the simplified case of a one-dimensional periodic 
lattice with spacing a, for which the RDF (3) will take the 
form 



1/2^,1 _ (r-naf 

> : e 2t7 2 n r > 0. 



a V it Jn 

n = l v 



Let us now calculate the Fourier transform of the function 
/i(r) = g(r)exp(— ra/cr 2 ): 



h{k) 



-n(q 2 +fc 2 cr 4 ) 
e 2tr 2 



a 2 +fc 2 cr 4 

e 2<t 2 - 1 



We can view h(fc) as a meromorphic function of complex 
variable fc. This function has poles at points k m = k' m + ik'^ e 
C which satisfy 

a 2 + k 2 ct 4 = 4 Turner 2 m e Z. 

Solving this equation and considering only poles with posi- 
tive imaginary parts, one gets 



ma I I 16m 2 Ti 2 a 4 
kL = — r V -1 + 1/1 + 3 and 



/2| 



V2a 2 



1+11 + 



16m 2 Ti 2 o" 4 



Now h(fc) can be presented as an inverse Fourier transform 
of h(k) and expanded into a series using the Cauchy residue 
theorem: 




2ni 



\k" cos(fc' r) - k' sin(fc' r)l 

|_m v m y m v m y J 



Since fc^ monotonically increases with m, in the limit r — » oo 
the term corresponding to m = 1 will be dominant, therefore 
asymptotically 

2 A e ra/a 2 -k"r 

S^) =- + [k'( cos^r) - k[ sHk[r)] 



a + a 2 \k,\ 2 



+ 



^ral<j 2 —k'lr 



)■ 



The exponent in the second term gives us the sought corre- 
lation length: 

-l 

r 2 



£7 



k'( - a/a 2 a 



- + l/-+4Tt 2 ^r -1 



Expanding this expression in powers of cr/a, one gets the 
following approximate formula: 

a 3 5cr 2 r c i 2L 

Z=-^ + — + a0\(cr/af\K — KL/5. 
2n A a A 2a L J u 
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